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The macroscopic quantum coherence in a biaxial antiferromagnetic molecular magnet in the pres- 
ence of magnetic field acting parallel to its hard anisotropy axis is studied within the two-sublattice 
model. On the basis of instanton technique in the spin-coherent-state path-integral representation, 
both the rigorous Wentzel-Kramers-Brillouin exponent and preexponential factor for the ground- 
state tunnel splitting are obtained. We find that the quantum fluctuations around the classical 
paths can not only induce a new quantum phase previously reported by Chiolero and Loss (Phys. 
Rev. Lett. 80, 169 (1998)) , but also have great influnence on the intensity of the ground-state 
tunnel splitting. Those features clearly have no analogue in the ferromagnetic molecular magnets. 
We suggest that they may be the universal behaviors in all antiferromagnetic molecular magnets. 
The analytical results are complemented by exact diagonalization calculation. 

PACS numbers:03.65.Bz, 75.45. +j, 75.50.Ee, 75.50.Xx 



I. INTRODUCTION 

In recent years, owing mainly to the rapid advances 
both in new technologies of miniaturization and in highly 
sensitive SQUID magnctomctry, there have been consid- 
erable theoretical and experimental studies carried out on 
the nanometer-scale magnets |l],|2| which have been iden- 
tified as candidates for the observation of macroscopic 
quantum phenomena (MQP) such as the tunnel- 

ing of the spin out of metastable potential minimum 
through the classically impenetrable barrier to a sta- 
ble one, i.e., macroscopic quantum tunneling (MQT), or, 
more strikingly, macroscopic quantum coherence (MQC), 
where the spin coherently oscillates between energeti- 
cally degenerate easy directions over many periods. In 
the semiclassical spin-coherent-state path-integral theory 
H , MQC is connected with the presence of a topological 
term in the Euclidean action Se{0,4>) arising from the 
nonorthogonality of spin coherent states, which is called 
Berry phase or the Wess-Zumino, Chcrn-Simons term: 
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the system, (9, </>) are polar and azimuthal spin angles, 
repectively. 

One of the manifestations of MQC is the ground-state 
tunnel splitting of magnetic systems. In the absence 
of an external magnetic field, it has been theoretically 
demonstrated that the ground-state tunnel splitting is 
completely suppressed to zero for the half-integer total 
spin ferromagnets or antiferromagnets with biaxial crys- 
tal symmetry 1^,0, resulting from the destructive inter- 
ference of the Berry phase in the Euclidean action be- 
tween the symmetry-related tunneling paths connecting 
two classically degenerate minima. Such destructive in- 
terference effect for half-integer spins is known as the 
topological quenching. But for the integer spins, the 
quantum interference between topologically different tun- 
neling paths is constructive, and therefore the ground- 



state tunnel splitting is nonzero. 

While such spin-parity effects are sometimes related to 
Kramers degeneracy, they typically go beyond this the- 
orem in rather unexpected ways. In the presence of an 
external magnetic field, as pointed out by Garg |^[| , the 
ground-state tunnel splitting can oscillate as a function of 
the field which is applied along its hard anisotropy axis in 
ferromagnets with biaxial crystal symmetry, and vanishes 
at certain values. This prediction is confirmed in one re- 
cent experiment carried out by Wernsdorfer and Sessoli 
Jl0| | . They developed a new technique to measure the 
very small tunnel splitting on the order of 10 _8 K in ferro- 
magnetic molecular Feg clusters. Indeed, they observed 
a clear oscillation of the tunnel splitting as a function 
of the magnetic field applied along the hard anisotropy 
axis, which is direct evidence of the role of the topo- 
logical spin phase (Berry phase) in the spin dynamics 
of these molecules. Although this field induced oscilla- 
tion's behaviour is investigated in great detail in ferro- 
magnetic systems now ]^-[l0|] , it is still less understood in 
antiferromagnetic molecular magnets such as Feio, Fe6, 
Vg, and antiferromagnetic ferritin ]ll]Jl^ . Golyshev and 
Popkov first studied MQC in a uniaxial antiferromag- 
netic fine particle in the presence of magnetic field [ pd| , 
and found similar oscillation behavior. But they only 
calculated the Wentzel-Kramers-Brillouin (WKB) expo- 
nent in weak field approximation, and paid no attention 
to its preexponential factor. Later, in 1998, Chiolero 
and Loss considered the oscillation's properties of a ring- 
like molecular magnet using an anisotropic nonliner a 
model (NLsM) [12[. In addition to the usual topolog- 
ical spin phase (Berry phase) term, they found a new 
quantum phase arising from fluctuations which is never 
seen in ferromagnetic molecular magnets. It is really a 
striking quantum property in antiferromagnetic molecu- 
lar magnets. Unfortunately, the fundamental physics of 
this novel quantum phase is less explored so far. 

In this paper, We would like to study the MQC of 
a biaxial symmetry antiferromagnetic molecular magnet 
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based on the two-sublattice model Jl3|~|T5t. By apply- 
ing the instanton technique in the spin-coherent-state 
path-integral representation Jl6| , we obtain the rigorous 
instanton solutions and calculate both the WKB expo- 
nent and preexponcntial factor in the ground-state tun- 
nel splitting. We will show that the quantum fluctua- 
tions around the classical paths can not only induce a 
new quantum phase previously reported by Chiolero and 
Loss 12 1, but also have great influence on the intensity of 
the ground-state tunnel splitting. Those features clearly 
have no analogue in the ferromagnetic molecular mag- 
nets. We suggest that they may be universal behaviors 
in all antifcrromagnetic molecular magnets. Due to the 
instanton methods are semiclassical in nature, i.e., valid 
in large spins and in continuum limit, we perform exact 
diagonalization calculations and find that they agree well 
with the analytical results. 



II. INSTANTON CALCULATIONS FOR 
SPIN-COHERENT-STATE PATH INTEGRALS 

We consider the ring-like antiferromagnetic molecular 
magnets(i.e. Feio, Fe6 and Vg) composed of N = 2n 
spins s regularly spaced on a circle lying in the xy-plane 
with an antiferromagnetic exchange interaction between 
them [ fH2|]l7f| . In general, the crystalline anisotropy at 
each site has biaxial symmetry. As usually for antiferro- 
magnets [^3|-^5|, we decompose the local spins into the 
two magnetic sublattices: S\ and S 2 with the same spin 
value S — ns. Then, the molecular magnet in an exter- 
nal magnetic field H acting along its hard anisotropy axis 
can be described by a spin Hamiltonian of the type 



U = jSx ■ S 2 + ( k 1 Sf z + k 2 Sf y - g^i B HS 1: 



k 2 S 2 2y 



gp B HS 2 



(1) 



where g is the lande factor, and /is is the Bohr magneton. 
k\ > k 2 > are the crystalline anisotropy coefficients, 
and we take the easy, medium, and hard axes as x, y, 
and z respectively for each sublattice. j is the exchange 
energy. In accordance with experimental results it will 
be assumed that j > k\ , k 2 for the strong antiferromag- 
netic coupling. Note that our two-sublattice configura- 
tion is only valid for the magnetic field H ^ H a . Here, 

H n = is the critical field at which the strong antifer- 

a 9Hb ^ b 

romagnetic exchange interaction jSi ■ S 2 is comparable 
to the Zeeman term g/iBH(Si z + S 2z ). 

In the semiclassical approach ]16| ], in order to obtain 
the ground-state tunnel splitting, one should compute 
the imaginary-time propagator in the spin-coherent-state 
representation: 



{hf \exp [-TLT]\fn) 
= J Vflexp(-S E ) 



= J V{0 1 }V{e 2 }V{4> 1 }V{4> 2 }exp(- J drC 



(2) 



over all trajectories which connect the initial state | hi) to 
the final state | hf). Here Oj, cf>j (j = 1,2) are the polar 
and azimuthal angles of each sublattice spin vector, the 
Lagrangian C include two parts pj: 



Co = ^ iS(l — cos6j)tj>j (t) 



(3) 



and 



C\ = J[l + cos#i cos 6 2 + sin 0i sin# 2 cos (</>i — (f> 2 )] 
+ (Kx cos 2 6j + K 2 sin 2 9j sin 2 fa) 

— g^sSH cos6 



j 3 , 



(4) 



3=1,2 



corresponding to the Berry phase term and the total Eu- 
clidean energy term E(8i,(j)i,9 2 ,(j) 2 ). Here, we have in- 
troduced K Y = hS 2 , K 2 = k 2 S 2 , and J = jS 2 . All 
terms in (^) are of apparent physical meaning. The 
first term is the exchange interaction energy, the sec- 
ond term is magnetic anisotropy energy and the third 
term is Zeeman energy. The dominant contribution to 
the imaginary-time propagator comes from finite action 
solutions of the Euler-Lagrange equations of motion (in- 
statons), which can be expressed as 



SS E 
SSe 



= 0, 



(5) 
(6) 
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(j = 1,2) denote the classical paths. 



where 

According to the instanton technique in the spin- 
coherent-state path-integral representation fl6[| , the in- 
stanton's contribution to the tunnel splitting A (not 
including the geometric phase factor generated by the 
Berry phase term in the Euclidean action) is given by 



A = 



p UJp 



S c l 

2^ 



1/2 



(7) 



where tu p is the small-angle precession or oscillation fre- 
quency in the well, and S c i is the classical action or the 
WKB exponent determined by Eqs. (||) and (Q) . Thepre- 
exponential factor po originates from the quantum fluctu- 
ations around the classical paths, which can be evaluated 
by expanding the Euclidean action to second order in the 
small fluctuations. 



1. Wentzel- Kramers- Brillouin exponent 

In our case, only low-energy trajectories with almost 
antiparallel S\ and S 2 contribute the path integral. It is 
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therefore, safe to say that tunneling of §2 follows tun- 
neling of Si |L4|]. For that reason we can replace #2 
and 4>2 by 7r — Oi + eg and ir + <f>i — e$ respectively 
(with \eg\ , \e^\ <C 1) in C. In the new coordinates, the 
imaginary-time propagator of the system can be repre- 
sented as 



V{eo}T>{e^ J T>{6i}V{ct>i} x 



exp 



drC(6\,(j)\,S0,e l i > ) 



(8) 



By simple algebra, up to the second order approximation 
about eg and e^, we obtain 

£ = i2S<p + 2 (Ki cos 2 9 + K 2 sin 2 9 sin 2 4>) 
-iS<j> sin 9 - Ki sin26> 

- K2 sin 29 sin 2 <j> — g[i B SH sin 9ee) 

- (K2 sin 2 9 sin 20) e^ 
-iS [(1 + cos 9) — sin Osg] e^ 

- (Agge^ + Ag^ege^ + A^e^) , 



(9) 



where 



S ■ J 

Agg = — i — COS 9(j) + — — K 1 cos 29 

g^sSH 



A 

A,/, 



-{-K2 cos 29 sin 2 1 
= i^ 2 sin26»sin2<?(>, 

= - sin 2 9 + K 2 sin 2 9 cos 2<p. 



cos( 



(10) 



In Eqs. (|9|) and (|10|), we have dropped the subscript of 9i 
and 4>i for clarity. Upon Gaussian integrating (|^) over eg 
and £0 one can obtain the following effective Lagrangian 



C ef f = i2S<p + 2 (Ki cos 2 9 + K 2 sin 2 9 sin 2 </>) + 



2J 



sm 



(11) 



Note that magnetic field enters only through the last 
term in Eq. (0) and has no influence on the tunneling 
barrier. Because of the condition Ki > K2, the equi- 
librium orientations of Si are (9,<p) = (f ,0) and (§,7r) 
which correspond to two degenerate classical minima of 
the energy, E = 0. It is obvious from symmetry that 
there are two different type instanton trajectories of op- 
posite windings around hard anisotropy axis. We denote 
them as ± instantons: 



<f> = 



±7r/2 



±7i7T (n = 1, 3, 5, . 



(12) 



To execute the first, we should seek the classical 
path(or paths) tt c i (r) = [9 (r) , 4> (r)) connecting the 
two minima, that minimizes the action Se = J drC e ff. 
This path satisfies the Euler-Lagrange equations of mo- 
tion (see also Eqs. (JsJ) and @) 



dC 



eff 



dC 



eff 



dQd (r) 



= 0. 



(13) 



Substituting the effective Lagrangian into Eq. (|l3|), we 
obtain 



d 



tPd9_ 

-2Ki + 2K 2 sin 2 $) sin 29 



. dS \ 
2J Shl2d [d^' l9 ^ H 



d_ 



J 



dr 



= 2K 2 sin 2 6'sin20. 



(14) 



(15) 



Consequently, a quasiclassical tunneling of Si may occur 
in xy-plane 9 = J, and then, Eq. ( p"5| ) reduces to sine- 
Gordon equation, 



or equivalently 



'dr 2 



dr 



io\ sin 2(f), 



■ loi sm ( 



(16) 



(17) 



1 /2 

where lui = ( 4 J ^ 2 ) . Under the boundary conditions 
in which the classical path approach the two minima as 
r — > ±00, we obtain an exact solution of this equation 



4> (t) = 2 arctan (exp (loit)) . 



(18) 



It is easily verified that (f> — > 0,7r, as r — * ±00. The 
corresponding classical action, i.e., the WKB exponent 
in the rate of quantum tunneling at finite magnetic field, 
can be evaluated by integrating the Euclidean action with 
above classical trajectories, and the result is found to be 



Re Sri ± i Im S r 



with 



ReS c , =45 



K2 
J 



ImSri = 2-kS 1- 



1/2 



H 



(19) 

(20) 
(21) 



where the positive and negative sign in Eq. ( |l9| ) are 
corresponding to ± instantons, respectivel y. 

It is clearly seen from Eqs. (|0|) and (|2l|), the clas- 
sical action has two unusual features in the presence of 
magnetic field. First of all, the real part of action has 
no dependence on the magnetic field and is determined 
by material parameters of the system only. This fea- 
ture is quite different from that in ferromagnetic molec- 
ular magnets, and can be understood easily from Eq. 
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(pi]), since the tunneling barrier remains unchanged un- 
der the magnetic field. Further, as shown from Eq. (21), 
if we ignore the contribution from quantum fluctuations 
around the classical paths, the ground-state tunnel split- 
ting which is proportional to exp(— Rc S c i) |cos(Im5 c j)| 
oscillates as the field H is increased, and the tunneling is 
thus quenched whenever 



(2g-n-l/2) 
H = 25 



(22) 



where n = 0,1,2,.... It is interesting to note that this 
result agrees well with Eq. (|ll]) in Ref. || found by Garg 
for ferromagnetic molecular Fes clusters if one makes the 
replacement J = 25 and sets A = 0. 



2. preexponential factor 

The second major step is to evaluate the preexponen- 
tial factor of small fluctuations around the classical in- 
stanton paths. We write 



9 (r) = 9 (r) + 59 (r) , 
0(r) =4>{t)+84>{t).. 



(23) 
(24) 



and evaluate the action to second order in (59, 5(f>). Writ- 
ing Se — S c i + 5 2 S, we have 



5 2 S 



dr 



21 

2J 



(gUsSHY 



U! 2 cos 2<j> ± i2g[i B SHuj 1 sin 4\ 59 2 



(u) 2 cos 



24) 54?), 



(25) 

1 /2 

where ujq = ( 4 J <^ 1 ) . Note that the 9 and 4 fluctua- 
tions are decoupled in Eq (p5|), and in the ^-fluctuation, 
an unusual term ±i2g^iBSHuji sin 4 distinguishing + in- 
stantons from — instantons appears. As we will show 
below, this extra term has important consequences at 
the high magnetic field. 

Now, the imaginary-time propagator is found to be 



(n f \ex P [-HT)\ hi) = exp (5±) D 



(26) 



with 



V{59} 



exp 



—oj 2 + lo\ cos 24 ± i2g[iBSHuj\ sin 4) 59 2 ~\ } 
D54, =^J x 

b 2 + (uj\ cos 24) 5q 



exp 



/ dT — 

2J 



(27) 



(28) 



where Ng and N<j> are the normalization factors. The 
fluctuation determinant for 4 is standard. Following the 
Eq. (2.44) in Ref. Q, we obtain 



5V 
Jtt 



1/2 



= 2wi 



/ Re5 c 
V 2?r 



1/2 



(29) 



For the ^-fluctuation determinant we find in the first or- 
der perturbation theory (The detailed calculation of Df g 
will be reported elsewhere.), for the high magnetic field, 



D 



so 



exp 



7 



1/2 



(30) 



where r] 



1 /2 

^SH = if (^f ) 18 USed aS a Sma11 P a " 

rameter in the high magnetic field. It is clearly shown 
in Eq. (p0|), the existence of magnetic field can bring 
a phase shift which approaches approximately ^ in the 
high-field regime. Note that the value of phase shift 
agrees well with that found by Chiolero and Loss Jl2| . 
On the other hand, as the field decreases down to zero, 
and thus r\ — > 00, the phase shift vanishes despite the 
breakdown of our first order perturbation calculations in 
the low-field regime. 

Combing the classical action and two fluctuation de- 
terminants, we arrive at the desired ground-state tunnel 
splitting, 



A = 4 exp 



4/2 



1/2 



/Re5 c 



V( 1 + f^ 
exp(-Re5 ci ) |cos$(#)l , 



V 2tt 



1/2 



where 



<!>.//; =2;,-5 ( i-JL 



¥ 2 v 



1/2- 



(31) 



(32) 



III. RESULTS AND DISCUSSIONS 

The semiclassical analysis presented so far applies 
strictly speaking only to a sizable number of spins with 
S> 1. However, as is often the case with such methods 
the results are valid (at least qualitatively) even down 
to a few spins of small size. This expectation is indeed 
confirmed by exact diagonalization calculations which we 
have performed on Hamiltonian (Q). Result for 5 = 5, 
and for some typical values ki = 0.03 K, = 0.01 K 
and j — 1.0 K is presented in Fig. 1, the critical field is 
found to be 7.44 T. Here, the units for the energy and 
magnetic field are taken to be Kelvin and Tesla, respec- 
tively. We can see that the numerical and semiclassical 
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approach show reasonable agreement in the whole mag- 
netic field regime. Since our perturbed calculation for 
the ^-fluctuation determinant is only valid in the high 
magnetic field, the agreement in the low-field regime is 
surprising in some ways. 

As shown in Fig.l, the ground-state tunnel splitting 
vanishs at the field ■§- = 1.0. This disappearance is ev- 
ident for the extra ? phase shift, since according to Eq. 
( f2l| ) , there should be a peak in usual. It is worthy noting 
that the extra ^ phase shift is not limited to our biax- 
ial symmetry antiferromagnetic molecular magnets case. 
Indeed, for the strong antiferromagnetic coupling, upon 
Gaussian integrating (||) over the small displacements eg 
and £(/,, one can obtain the effective Lagrangian in general 
form 



C e ff — i2S(f> 



-El 



52 
2J 



, (33) 



where the detailed form of E (9, (j>) depends on the 
system investigated. Then, an extra term similar to 
H2g/j,BSHu!ism(j) in Eq. ( p5|) will appear after we ex- 
pand the small fluctuations around the classical instanton 
paths to the second order, and gives the ? phase shift. 
Thus, we conclude that the \ phase shift induced by 9- 
fluctuation may be a universal quantum behavior in all 
antiferromagnetic molecular magnets systems. 

In the figure, another interesting feature is the peak 
height of the ground-state tunnel splitting. As the mag- 
netic field increases, the peak height first drops signifi- 
cantly in the low-field regime, and then keeps invariant 
up to the critical field H a . This may be qualitatively 
understood from Eq. (pl|). Because the WKB expo- 
nent has no relevance with the magnetic field, the most 
essential dependence of the ground-state tunnel split- 
ting on the field comes from the preexponential factor 



Po 



4exp 



1+ 



which undergoes a dramati- 



cs • 



cally change only in the low-field regime. 

At the end of this section, in order to support the ex- 
perimental relevance of our results, we give some esti- 
mates for the ferric wheel, Feio, f° r which N — 2n = 10, 



4T,a = 0.03 



s — I, S — ns — 12.5. If one takes - 2 — 
and y- = 0.01 as the typical parameters values |lj], then 
the simple algebra demonstrates that the ground-state 
tunnel splitting has 2S = 25 oscillations of magnitude 
A k 1.3 K and period 4 T. From Eq. (|o|), the phase 
shift will be visible when ^77 ~ 3, or H = 10 T. There- 
fore, all quantities appear to be well within experimental 
reach. 



IV. CONCLUSION 

In summary, we have investigated the MQC phe- 
nomena in biaxial symmetry antiferromagnetic molec- 



ular magnets. Our discussion is based on the two- 
sublattice model that includes anisotropy and magnetic 
field. On the basis of instanton technique in the spin- 
coherent-state path-integral representation, both the rig- 
orous Wentzel-Kramers-Brillouin exponent and preexpo- 
nential factor for the ground-state tunnel splitting are 
obtained. We have outlined here two prominent features 
in our tunneling scenario: (?) In addition to the usual 
topological term in the classical action, a new quantum 
phase arising from the quantum fluctuations around the 
classical paths is found to contribute to the tunneling 
oscillations. This result coincides with that reported by 
Chiolero and Loss p2| . (ii) The magnetic field appears 
to have no influence on the tunneling barrier. Thus the 
main dependence of the tunneling peak height on the 
field comes from the quantum fluctuations, this leads to 
a sudden drop of peak's height in the low-field regime. 
Both two features clearly have no analogue in the fer- 
romagnetic systems |§H1QI- We suggest that they may 
be universal behaviors in all antiferromagnetic molecular 
magnets. 

We realize that our result is based on the instanton 
method which is semiclassical in nature, i.e., valid in 
large spins and in continuum limit. We perform exact 
diagonalization calculation to check its validity, and find 
that it agrees well with the analytical result in the regime 
where a comparison is possible. 

Recent experiments have rekindled interest in the field 
of quantum tunneling of the molecular magnets.. Most 
notable has been the discovery of resonant quantum tun- 
neling between spin states in the ferromagnetic system of 
spin-10 molecules such as Mn^Ac jy and Fes 0- Since 
the antiferromagnetic molecular magnets are proposed as 
better candidates for observing the phenomena of MQT 
and MQC |D| compared with the ferromagnetic ones, we 
hope that our predictions on antiferromagnetic molecular 
magnets can be confirmed in experiments in the future. 

Acknowledgments 

The financial supports from NSF-China (Grant No. 
19974019) and China's "973" program are gratefully ac- 
knowledged. 



[1] R. Sessoli, D. Gatteschi, A. Caneschi, and M. A. Novak, 
Nature 365, 141 (1993). 

D. Gatteschi, A. Caneschi, L. Pardi, and R. Sessoli, Sci- 
ence 265, 1054 (1994). 

L. Thomas, F. Lionti, R. Ballou, D. Gatteschi, R. Sessoli, 
and B. Barbara, Nature 383, 145 (1996). 
J. M. Hernandez, X. X. Zhang, F. Luis, J. Tejada, J. R. 
Friedman, M. P. Sarachik, and R. Ziolo, Phys. Rev. B 
55, 5858 (1997). 

D. A. Garanin and E. M. Chudnovsky, Phys. Rev. B 56, 
11102(1997). 



5 



[2] A. -L. Barra, P. Debrunner, A. Catteschi, C. E. Schulz, 
and R. Sessoli, Europhys. Lett. 35, 133 (1996). 
C. Sangregorio, T. Ohm, C. Paulsen, R. Sessoli, and D. 
Gatteschi, Phys. Rev. Lett. 78, 4645 (1997). 

[3] For an overview, see e.g., Quantum Tunneling of Magne- 
tization - QTM'94, edited by L. Gunther and B. Barbara 
(Kluwer, Dordrecht, the Netherlands, 1995), and refer- 
ences therein. 

[4] E. M. Chudnovsky and J. Tejada, Macroscopic Quan- 
tum Tunneling of the Magnetic Moment, Cambridge Un- 
verisity Press, 1998. 

[5] E. Fradkin, Field Theories of Condensed Matter Systems, 
Addison- Wesley, New York, 1991. 

[6] D. Loss, D. P. Divincenzo, and G. Grinstein, Phys. Rev. 
Lett. 69, 3232(1992). 

[7] J. Von Delft and C. L. Henley, Phys. Rev. Lett. 69, 
3236(1992). 

[8] A. Garg, Europhys. Lett. 22, 205(1993). 

[9] A. Garg, Phys. Rev. B 60, 6705(1999). 
[10] W. Wernsdorfer and R. Sessoli, Science 284, 133 (1999). 
[11] V. Yu. Golyshev and A. F. Oppkov, Europhys. Lett. 29, 
327(1995). 

[12] A. Chiolero and D. Loss, Phys. Rev. Lett. 80, 169(1998). 
[13] B. Barbara and E. M. Chudnovsky, Phys. Lett. A 145, 
205(1990). 

[14] E. M. Chudnovsky, J. Magn. Magn. Mater. 140-144, 
1821(1995). 

[15] H. Simanjuntak, J. Phys. Cond. Matter 6, 2925 (1994). 
[16] A. Garg and G. -H. Ki m, Phys. Rev. B 45 , 12921(1992). 



[17] A. K. Zvezdin, eprint, cond-mat/0004074 , (2000) 



Figure caption 

Fig.l. ^ versus with S = 5, k x = 0.03 K, 
k 2 = 0.01 K, j = 1.0 K and H a = 7.44 T. Here, the units 
for the energy and magnetic field are taken to be Kelvin 
and Tesla. The solid line and symbols represent, respec- 
tively, the ground-state tunnel splitting predicted by the 
semiclassical instantons approach and that obtained by 
the exac diagonalization method. 
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Fig.l 

Macroscopic quantum coherence 
H. Hu et al. 




